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ABSTRACT 



In this paper, the compiled long-term optical and infrared measurements of some 
blazars are used to analysis the variation properties and the optical data are used to 
search for periodicity evidence in the lightcurve by means of the Jurkevich technique 
and the discrete correlation function (DCF) method. Some explanations have been 
discussed. 

Subject headings: Active Galactic Nuclei (AGNs): variability - periodicity 
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1. Introduction 

The nature of the central regions of quasars and other Active Galactic Nuclei (AGNs) is still 
an open problem. The study of AGN optical variability can yield valuable information about the 
mechanisms operating in these sources, with important implications for quasar modeling (see, 
for instance. Fan et al. 1998). Some particular objects have been claimed to display periodicity 
in their lightcurves over a variety of timescales (e.g. Jurkevich 1971; Hargen-Thorn et al. 1987; 
Sillanpaa et al. 1988; Webb et al. 1988; Babadzhanyants & Belokon 1991; Kidger et al. 1992; Liu 
et al. 1995; Marchenko et al. 1996; Fan et al. 1997, 1998, Fan & Su 1998; Fan 1999a; Zhang et 
al. 1998; Fan & Lin 1999a; Fan et al. 1999a; Lainela et al. 1999; Lin 1999; Su 2000 and references 
therein), but, in general, the clear identification of periodic behaviour has been very elusive due to 
the complexity of the optical lightcurves and the lack of databases large enough as to provide an 
adequate sampling over large periods. 

Blazars are an extremely subclass of AGNs, characterized as rapid and high amplitude 
variations, high and variable polarization, superluminal radio components, nonthermal continuum. 
Blazars have been monitored for a long time, compilations are available for optical (Fan & Lin 
1999b, c) and infrared (Fan & Lin 1999d; Fan 1999b) bands. The infared data are available in 
http: / /xxx.lanl.gov/astro-ph/990810^ and |http:// xxx.lanl.gov/astro-ph/9910269|. The long-term 



measurements make it possible for one to search for periodicities from the light curves and 
to discuss the long-term variability properties. Because the characteristic of the astronomical 
measurements, we adopted the Jurkevich (Jerkevich 1971) and DCF (Discrete Correlation 
Function) (see Fan & Lin 1999a) methods to the available data in searching for the periods. In 
the second section we will give the variation properties, in the third section, we present the results 
and finally in section 4, we give some discussions. 
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2. Variation Property 

Prom the complication (Fan et al. 1999b,c,d; Fan 1999b), we found that the largest variations 
are comparable in both the optical and infrared bands. Largest variations at different wavelengths 
increase with decreasing wavelength. The variations are also correlated with the highest observed 
optical polarization, which can be explained using the beaming model. Some objects show that 
the spectra flatten when the sources brighten while the opposite cases are for some others. For the 
largest infrared and optical variations, we listed them in Table 1. 

3. Periodicity Analysis Methods 

3.1. Jurkevich Method 

The Jurkevich method ( Jurkevich 1971, a Iso see Kidger et al. 1992; Fan et al. 1998a; Fan 
1999a) is based on the expected mean square deviation and it is less inclined to generate spurious 
periodicity than the Fourier analysis used by other authors (e.g. Fan et al. 1997). It tests a run 
of trial periods around which the data are folded. All data are assigned to m groups according 
to their phases around each trial period. The variance for each group and the sum of all 
groups are then computed. If a trial period equals the true one, then reaches its minimum. 
So, a "good" period will give a much reduced variance relative to those given by other false trial 
periods and with almost constant values. A further test is the relationship between the depth of 
the minimum and the noise in the "fiat" section of the curve close to the adopted period. If 
the absolute value of the relative change of the minimum to the "flat" section is large enough as 
compared with the standard error of this "flat" section (say, five times) , the periodicity in the 
data can be considered as significant and the minimum as highly reliable (Fan et al. 1998a). 
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3.2. DCF Method 

The DCF (Discrete Correlation Function) method is intended for analyses of the correlation 
of two data set. It is described in detail by Edelson & Krolik (1988) (also see Fan et al. 1998b). 
This method can indicate the correlation of two variable temporal series with a time lag, and can 
be applied to the periodicity analysis of a miiquc temporal data set. If there is a period, P, in the 
lightcurve, then the DCF should show clearly whether the data set is correlated with itself with 
time lags of r = and t = P (see Fan k. Lin 1999a). We have implemented the method as follows. 

Firstly, wc have calculated the set of unbinned correlation (UDCF) between data points in 
the two data streams a and 6, i.e. 

where Oj and hj are points in the data sets, a and b are the average values of the data sets, and (Ta 
and cTft are the corresponding standard deviations. Secondly, we have averaged the points sharing 
the same time lag by binning the U DCFy in suitably sized time-bins in order to get the DCF for 
each time lag r: 



DCF{T) = ^^UDCFijiT), (2) 



where M is the total number of pairs. The standard error for each bin is 



^^^^ = W^^^ ^^^'^^'^ ~ Z^CF(r)]2}0 ^ (3) 



Results and Discussion 



Very recently, we compiled four data bases for BL Lac objects and OVV/HPQs in the optical 
and infrared bands (Fan & Lin 1999b,c,d; Fan 1999b). The data bases can be used to discuss the 



-6- 



long-term variability properties, the correlated variations and to search for periodicity signatures 
in the light curves. When the both methods are adopted to those optical data, the both methods 
give consistent results. The results are presented in Table 2. It is interesting to notice that there 
is a common period of ~ 1.0 year, which is likely from the effect of the Sun on the measurements. 
The detail analysis is presented in the paper by Fan et al. (1999c) 

For the long-term periodicity of variations, there are several explanations: the double black 
hole model, the hectic jet model, the slim disk model, and the effect of external perturbations to 
the accretion disk. (e.g. Sillanpaa et al. 1988; Meyer & Meyer-Hofmeister 1984; Horiuchi & Kato 
1990; Abraham & Romero 1999; Fan et al. 1997, 1998a, 1999; Villata & Raiteri 1999; Romero et 
al. 2000). 

This work has been supported by the National PanDeng Project of China and the National 
Natural Scientific Foundation of China. I thank Prof. K.S. Cheng to provide me with the 
opportunity to present this work and Prof. R.G. Lin and C.Y. Su for their help. 
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Table 1: The largest variations of Blazars 



Name 


Airiopt 


AmiR 


Popti%) 


Name 


Airiopt 


AmiR 


Popti%) 


(1) 


(2) 


(3) 


(4) 


(1) 


(2) 


(3) 


(4) 


0048-097 


2.7 


6.55 


27.2 


0109+244 


3.07 


1.58 


17.3 


0118-272 


1.05 


0.67 


17. 


0138-097 


1.52 


1.69 


29.3 


0215+015 


5.0 


2.69 


20. 


0219+428 


2.0 


1.61 


18.0 


0235+164 


5.3 


5.0 


44. 


0323+022 


1.3 


0.97 


10.4 


0420-014 


2.8 


2.88 


20. 


0422+004 


2.2 


3.25 


22. 


0521-365 


1.4 


1.25 


11. 


0537-441 


5.4 


3.0 


18.8 


0716+714 


5.0 




29.0 


0735+178 


4.6 


2.47 


36. 


0736+017 


1.35 


2.71 


6. 


0754+100 


3.16 


1.88 


26. 


0818-128 


3.78 


2.14 


36. 


0823-223 


1.41 


2.32 


11. 


0829+046 


3.58 


2.15 


12. 


0851+202 


6.0 


3.87 


37.2 


0912+297 


2.25 


2.27 


13. 


1101+384 


4.6 


2.03 


16.0 


1144-379 


1.92 


3.46 


8.5 


1147+245 


1.0 


1.34 


13.0 


1156+295 


5.0 


4.47 


28. 


1215+303 


3.3 


1.05 


14. 


1219+285 


3.13 


2.46 


10. 


1253-055 


6.7 


4.57 


44. 


1308+326 


4.17 


5.55 


28. 


1418+546 


4.8 


1.59 


24. 


1510-089 


5.4 


1.25 


14. 


1514-241 


3.0 


2.64 


8.0 


1538+149 


3.7 


0.89 


32.8 


1641+395 


3.0 


3.16 


35. 


1652+398 


1.30 


2.37 


7.0 


1727+502 


2.1 


0.82 


6.0 


1749+096 


2.7 


2.21 


32. 


1807+698 


2.0 


1.02 


12. 


1921-293 


2.64 


3.06 


17. 


2005-489 


0.53 


0.31 


2.0 


2155-304 


1.85 


1.88 


14.2 


2200+420 


5.31 


2.93 


23. 


2223-052 


5.0 


3.96 


17.3 


2240-260 




1.66 


15.1 


2251+158 


2.5 


1.57 


19. 


2254+074 


3.27 


1.36 


21. 



Notes to the table: 
Col. 1 The name of the source; 
Col. 2, The largest optical variation; 
Col. 3 The largest infrared variation; 
Col. 4, The largest optical polarization 
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Table 2: Investigation Results of Periodicity of 14 Blazars 



Name 


Pi (/)(year) 


P2 (/)(year) 


P3 (/)(year) 


P4 (/)(year) 


0219+428 


0.99 


2.14 


2.93 


4.52 


0235+164 


1.01 


1.56 


2.95+0.25 


5.87+1.3 


0735+178 


1.01 


14.39+0.51 


18.58+0.52 


28.62+0.58 


0754+100 


1.00 


17.85+1.3 


24.7+0.7 




0851+202 


1.01 


5.53 


11.75+0.5 




1215+303 


0.99 


2.05+0.1 


4.45+0.05 


6.89+0.34 


1219+285 


0.99 


1.97 


9.0+0.1 


14.84+1.55 


1226+023 


0.98 


2.00 


13.5+0.2 


22.5+.02 


1253-055 


1.2 


7.1+0.4 






1308+326 


1.00 


2.03 


2.95 


6.07+0.1 


1418+546 


0.99 


2.99 






1514-241 


1.00 


2.01 


4.09 


9.35+0.1 


1807+698 


0.99 


6.0 


18.18 




2155-304 




4.16 


7.0 




2200+420 


0.99 


13.7+0.4 


19.95+0.70 


29.6+0.30 



Notes to the Table: 

Pi, P2, P3, P4 are the possible periods derived from the light curves 
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